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Both the stretch and bend vibrational spectrum and the intensity of infrared tran- 
sitions in a tetrahedral molecule are studied in a U(2) algebraic model, where the 
spurious states in the model Hamiltonian and the wavefunctions are exactly removed. 
As an example, we apply the model to silicon tetrafluoride SiF4. 

I. INTRODUCTION 

In recent years, algebraic methods [1,2] have been developed to be an effective theory for descriptions 
of vibrations, rotations, and rotation- vibration interactions in polyatomic molecules [3,4]. U(4) and 
U(2) algebraic model have mostly been used so far in the analysis of experimental data. U(4) model 
[5,6] takes the rotation and the vibration into account simultaneously but becomes quite complex 
when the number of atoms in the molecules becomes larger than 4, while U(2) model is particularly 
well suited for dealing with the stretching vibrations of polyatomic molecules such as the octahedral 
and benzene-like systems [7,8]. Those two models are still extensively used for small molecules [9,10]. 
Recently, an active step in the development of algebraic models has been taken by Iachello and Oss 
to incorporate the bending modes into the models for benzene [11] and acetylene [12]. Extending 
U(2) algebraic model [3] to describe both the stretching and the bending vibrations in X3 [13] and X4 
molecules [14], Frank and his co-workers have introduced a symmetry-adapted algebraic model, and 
established an explicit connection between algebraic and configuration space calculations. Developing 
the local mode model [15], we have recently proposed an algebraic model of boson realization [16- 
18] for the complete vibrations by taking full advantages of the discrete symmetries of molecular 
systems. It is worth mentioning that U(k-l) and U(n) algebraic approach have been introduced for 
the k=3m-3 rotational and vibrational degrees of freedom of m-atomic molecules [19], and for n-1 
stretching vibrational degrees of freedom of polyatomic molecules [20], respectively. However, those 
two approaches are less feasible than U(4) and U(2) model in application. 

In order to study vibrational modes of medium-size molecules, one usually chooses a tetrahedral 
molecule as one of good samples. Its excited stretching vibrational states were explained in the local 
mode model [21]. In our previous papers, the vibrational spectrum of methane was analyzed in terms 
of bosonic operators [16] and q-deformed harmonic oscillators [23], however its infrared intensity was 
not taken into account. Recently new algebraic SU(2) model [24] has been proposed for tetrahedral 
molecules. This model needs at least eight parameters for the calculation of transition intensities. The 



lack of experimental data did not allow such a fit for the methane molecule, and interactions between 
the stretch and the bend were neglected in the calculation of energy levels. Leroy et al. have up to date 
presented an algebraic model [22] only for highly excited stretching modes and infrared intensities in 
this system. Although a precise infrared transition model has been constructed in algebraic scheme, 
it will be complex when the bending vibrations are considered. 

In this letter, both stretching and bending vibrations and intensities of infrared transition of all 
active modes in a tetrahedral molecule are studied in a U(2) algebraic model based on our new method 
for eliminating the spurious states and spurious components in Hamiltonian. Our model transition 
operator contains only five parameters. One of its applications is presented to silicon tetrafluoride 
SiF 4 . Calculated results demonstrate that our method for dealing with spurious states are also effective 
for U(2) algebraic model, which open the possibility to use this model for analyzing the complete 
vibrations in large molecules where spurious states exist. In Sec. II, considering the symmetry of a 
tetrahedral molecule, we set up the model Hamiltonian, in which the interactions related to spurious 
states are exactly removed. Its infrared intensity is studied in Sec. III. Application to SiF4 is given 
in Sec. IV, where symmetrized bases are used to simplify calculations. Concluding remarks are made 
in Sec. V. 

II. HAMILTONIAN 

For a tetrahedral molecule XY4, we introduce ten U(2) algebras to describe ten interactions between 
atoms: U/(2) (1 < j < 4) for X-Y and U M (2) (5 < \i < 10) for Y-Y interactions. The assignments of 
the Cartesian coordinate system are the same as in Ref. [16]. Each U Q (2) (1 < a < 10) is generated 
by the operators {N a , J+, a J-, a , Jo,a}, satisfying the following commutation relations: 

[Jo,a, J±A — ±^Q/3>/±,a, [J+,a, J-,p] = 2<5 Q /3 Jo,a, 

[N a ,J ,p}=0, [N a ,J±,p]=0. 
where N a is related with the Casimir operator of U(2): 

2io,a + J+,aJ-,a + j-,*J+,a = N a (N a /2 + 1). 

Denote by v a the number of quanta in the ath bond. The local basis states for each bond are labeled 
by the eigenvalue N a of N a and v a , and written as \N a , v a ). Their products provide the local bases: 

|JVi,«l)|iV2,«2)---|JViO,«10> = \N a ,v a ). 

Those Nj of equivalent bonds are equal to each other: Nj — N s , and — Nb, where and hereafter, 
the indexes j, fi, and a run from 1 to 4, 5 to 10, and 1 to 10, respectively, and footnotes s and b refer 
to the X-Y couplings and the Y-Y couplings, respectively. 

There are three kinds of 0(2) invariant combinations of those generators: 

Ha = (J+,uJ-,a + J-,a, J+,a)/2 — N a /2, 

H a ,0 = 2J ,Jo,p - N a N p /2, a^p, (2.1) 

V a ,0 = J+, a J-,/3 + J-,a J+,0, Cfj^ f3. 
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Their matrix elements in the local bases are given in Ref. [13]. The operators H^f correspond to the 
energy of the ath Morse oscillator. The operators # Q ,/3 describe anharmonic terms of the type v a Vfj, 
while the operators V a ^ describe interbond couplings which, in configuration space, are of the type 
r Q r^, where r a and are the displacement vectors of bonds a and (3 from their equilibrium positions. 

Considering that ^ J +iM (or J2J-^) relates only to spurious states [23], we obtain the following 
Td invariant Hamiltonian in terms of those three kinds of operators: 

H = H s + Hb + Hgi,, 

4 4 4 



H s — \ s \ ^ Hj 1 + \ S 2 ^ #,j + \ S 3 ^ Vij, 
10 10 

#6 = A b i Hff + ^62 ^ H^ v 

5 v— p,^Z,ii<v— 6 

7 7 

+ a 63 E K 

^l— 5 /i— 5 

H s b = ^sbl I ^2 (^.M _ ^l.M+s) + ( #2,5 - E #2,,u + #2,9 + #2,10 j 
1/4=5 \ /i=6 / 

+ ^2 (^3,2/j - #3,2/j-l) + ( --#4,5 - #4,6 + ^ #4,/i — #4,10 j \ 

+ A sb2 J ^ (y^ - V W ) + f^s-^^ + ^g + ^oj 

5 / 9 \ ^ 

+ Yl - ^3,2 M -l) + -14,5 - ft,6 + 2 ^ - 



(2.2) 



,10 

fi— 3 \ M— 7 

where # s , #{,, and £f s (, describe the stretching interaction, the bending one, and the interaction 
between the stretching and bending modes, respectively. The Hamiltonian preserves the quantum 
number V — v a ■ 



III. INTENSITIES OF INFRARED TRANSITION 



The infrared active mode is F 2 . The absolute absorption intensities from state v' to v are given by 

Pw = \{v\f x \v')\^ + \{v\f v \v')\^ + \{v\f z \v')\\ 

where v vv * is the frequency of the observed transition, T x , T y , and T z correspond to the three com- 
ponents of the infrared transition operator T, and the state \v) denotes \N a ,v a ) for short. All other 
constants are absorbed in the normalization of the operator T. The three components of T are 

f x = Cs (ti - i 2 + h -U) + Cb (4 -k) + Csb (ti + i 2 + h + U){k -ig), 
% = (s (ti - h - t 3 + U) + Cb {i 7 - iio) + Csb (h + h + h + h)(t 7 - i w ), (3-2) 
T z = Cs (ti + h - t 3 -U) + Cb (4 -h) + Csb (ti + h + h + U)(h -t s ), 
where Cs, Cb, an d Csb are parameters, and t a is the local operator for the crth bond. The term with 
Csb is one of the higher order contributions of T, which is necessary for describing both the stretching 
and the bending active modes. The matrix elements of t a are taken to be [11] 
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(N a ,v a \i a \N a ,v' a ) = exp(~i] a \v a - v' a \). (3.3) 

Those rjj for equivalent bonds are equal to each other: r/j = r] s , and i]^ = r] b . 

The five parameters in the transition operator of Eqs. (3.2)-(3.3) will be determined by fitting 
observed data. The calculated intensity can be used to check assignments and in the study of in- 
tramolecular energy relaxation in tetrahedral molecules. 

IV. APPLICATION 

We now apply this model to study the complete vibrational spectrum and the infrared intensity 
of SiF4. To our knowledge, there are sixteen observed vibrational energy levels and thirteen infrared 
intensity data for SiF4 [25]. 

At first, we calculate the Hamiltonian matrix. The calculation for energy levels will be greatly 
simplified if the symmetrized bases are used. For the stretch and bend states of tetrahedral molecules, 
the symmetrized bases were given in Ref. [16], where the spurious states in the bases were eliminated. 
For stretching overtones of large molecules, the symmetrized bases have been recently presented by 
Chen et al. [26,27]. In those symmetrized bases the Hamiltonian matrix is a block one. According to 
Morse potentials for the stretch and the bend vibrations in the boson-realization model [28] for this 
molecule, we take two boson numbers N s and N b to be 100 and 15, respectively. The parameters in 
Hamiltonian are determined by fitting the observed data, and given in cm -1 as follows: 

A s i = -3.311, A s2 = -4.331, A s3 = -0.573, A 6 i = 97.982, A b2 = -1.197, 
A b3 = 74.438, A b4 = -4.173, X sbl = -0.145, \ sb2 = 0.199. 

The experimental data and the calculated values are listed in Table I. The standard deviation in this fit 
is 1.188 cm -1 . From those parameters we can calculate the other values. It is worth mentioning that 
McDowell et al. [25] described the same energy levels by Dunham expansion with more parameters. 
This method does not provide explicitly wave functions so that some physical properties such as 
transition intensities arc hard to be calculated. 

TablcT 

Then, we compute the infrared intensity. Due to the Wigner-Eckart theorem, it is sufficient to 
calculate only the z component of the transition operator, T z , in the symmetrized bases. Fitting 
the thirteen observed infrared intensities, we determine the parameters in the transition operator as 
follows: 

r) a = 40.158, r) b = 22.853, ( s = 3.614, ( b = 3.286, ( sb = 64.959. 

In Table I we only list those calculated intensities to compare with known observed data. The other 
calculated intensities and energy levels can be obtained from us upon request. 

One may observe in Table I that most of the calculated intensities are in good agreement with the 
experimental values, but a few are not. Those differences may come from that the observed intensities 
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were only approximately accurate [25], and that the other higher order contributions to the operator 
T were neglected. It is worth pointing out that there are also differences in magnitude between 
the calculated and the observed intensities for the stretching vibrations of octahedral molecules in 
U(2) model [29] because the same simple matrix elements Eqs.(3.3) were used. The more accurate 
experimental data are needed to improve the model. 

V. CONCLUDING REMARKS 

We have used a U(2) algebraic model for studying the stretching and bending vibrations and infrared 
intensity of a tetrahedral molecule. The model Hamiltonian with nine parameters and the model 
transition operator with five parameters provide quite good fits to the published experimental data 
of silicon tetrafluoride with standard deviations 1.188 cm -1 and 1.775, respectively. This is based on 
our new methods for constructing symmetrized bases [16] and for removing both the spurious states 
in the wavefunction space and the spurious components in the Hamiltonian [23]. For comparison, we 
also studied this molecule in the boson-realization model, and obtained the corresponding standard 
deviations 1.985 cm -1 and 1.745 [28]. Through this example of applications of the model, we believe 
that our method for treating spurious states is useful for the model for other polyatomic molecules. 

It is shown that anharmonic resonances, such as Darling-Dennison and Fermi resonances, are very 
important in descriptions of highly excited vibrational states in molecules. Darling-Dennison reso- 
nances can be included in the model by adding higher-order terms of the operators of Eq.(2.1). Fermi 
resonances can be taken into account using perturbation theory. It should be pointed out that Fermi 
resonances can be easily included in the extended local mode model [30,31] and the boson-realization 
model [17], while they are described by the nondiagonal matrix elements of Majorana operators in 
U(4) algebraic model. The coupling parameters in Hamiltonian we have obtained can be related to 
force field constants to be used in conjunction with a kinetic energy operator in a Schrodinger equation 
[32]. We will consider the physical meaning of the parameters in future investigations. 
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Table I. Observed and calculated energy levels and relative intensity 



Obs.[25] 


Calc. 


r 


E (cm 1 


) Intensity 


V 


E (cm 1 ) Intensity 


E 


264.2 




1 


264.415 




F 2 


388.4448 


500 


1 


388.858 


500.059 


F 2 


776.3 


0.9 


2 


775.327 


2xl0~ 6 


M 


800.8 




1 


799.770 




F 2 


1031.3965 


5000 


1 


1029.677 


4999.867 


E 


1064.2 




2 


1064.186 




F 2 


1164.2 


1.4 


3 


1164.169 


3xl0~ 5 


F 2 


1189.7 


40 


2 


1188.631 


40.065 


F 2 


1294.05 


2.4 


2 


1293.903 


4 x 10~ 4 


F 2 


1418.75 


0.1 


2 


1418.533 


0.005 


F 2 


1804.5 


0.7 


3 


1804.706 


5xl0~ 5 


F 2 


1828.17 


7 


2 


1828.745 


3.888 


F 2 


2059.1 


1.2 


2 


2058.010 


3.744 


F 2 


2602.55 


0.007 


4 


2603.789 


4 x 10~ 8 


F 2 


2623.8 


0.015 


3 


2623.678 


0.001 


F 2 


3068.5 


0.015 


3 


3069.126 


lxlO" 5 
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